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Abstract 

We formally derive the hydrodynamic limit of a system modelling a 
bosons gas having a condensed part, made of a quantum kinetic and a 
Gross-Pitaevskii equation. The limit model, which is a two-fluids Euler 
system, is approximated by an isentropic system, which is then studied. 
We find in particular some conditions for the hyperbolicity, and we study 
the weak solutions. A numerical example is given at the end. 

1 Introduction 

There has been a lot of research done about Bosc-Einstcin condensates in math- 
ematics and physics, especially since the Nobel Prize 2001 was awarded to Carl, 
Cornell and Ketterle who succeeded in creating one for the first time. In par- 
ticular, mathematicians have been working on various models to see if the pre- 
dictions of the physicists could be justified mathematically. 

There are actually several types of models to describe a Bose-Einstein con- 
densate [13]. The only one which is valid without any approximation is the 
atomistic one, that is, the linear Schrodinger equation for the P bodies density: 

ihdt-^P = Hp-^P (1) 

where Hp is the Hamiltonian, i.e. the energy operator, of the system, and h is 
Planck's constant. Other models can be found depending on the regime under 
consideration. At zero temperature, the gas is entirely condensated. Assuming 
that it is in a non-dissipative trap, the one-body density is governed by the 
Gross-Pitaevskii equation: 

where U is the trapping potential. When the temperature of the gas is close to 
zero, the gas is composed of a condensate and a normal component. The conden- 
sate is still expected to be governed by a Gross-Pitaevskii-type equation, with 
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Figure 1: Various models deseribing a Bose gas, and their interactions. 

coupling terms taking into account the mass and momentum exchanges with 
the normal component. The latter can be described in a probabilistic way using 
kinetic equations; in this case, a quantum Boltzmann equation, first proposed 
by Nordheim [12], with additional terms. We could also use a fluid description 
for the normal part together with the usual Gross-Pitaevskii equation for the 
condensate; still another possibility is to use a fluid description for both parts 
of the gas, leading to a two fluids model as predicted by Landau [10]. The way 
these models are related is represented in Figure 1. 

In [14], a kinetic model for the non-condensate part has been derived from 
(1) using the BBGKY expansion in the low density limit. The condensated part 
is then described by a Gross-Pitaevskii equation, which takes into account both 
the interactions between the atoms of the condensate, and between them and 
those of the non-condensate part (called normal fluid). 

What we are interested in is to link this model with models at a larger scale, 
namely, fluid models. This can be done considering the fast relaxation limit 
for the kinetic equation, together with the semiclassical limit for the Gross- 
Pitaevskii equation. 

Our starting point is the following model: let 



ifi = ip{t, X, f ) > 



denote the phase-space density function of the normal fluid, and 



(2) 
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denote the wave function of the condensated (or superfluid) phase, with ps being 
the mass of the superfluid. Here, t G K+ is the time variable, and x, w G K.^ for 
A'' > 1 are the position and velocity variables. 

In [14], Pomeau, Brachet, Metens and Rica derive the following coupled 
system of equations: 

dt^p + I'.Va^v? - ^xVi.S/^ip = Q{ip) + psQi{ip) 

ih (3) 

ihdt'H = - yA^r* + V^2* + yQ2(V5)*- 

The first equation expresses a balance between the left-hand term, which is a 
Vlasov operator for a gas in the external field V^^T^, and the right-hand term, 
which comes from the normal fluid-normal fluid and normal fluid-supcrfluid 
interactions. The second equation is a cubic Schrodinger equation (because 
of the form (8) of V2), with an interaction term with the normal fluid. The 
first collision operator Q{ip) is the usual Boltzmann-Nordhcim operator for the 
bosons, modelling the collisions between two particles from the normal fluid. It 
writes 

Q{^)= I I (V'V:(1 + ^)(1 + ^')- W*(l + ^')(l + V':))di;*dn (4) 

with the usual notations 

if ^ ip{t,x,v), ip^ = lfi{t,X,V*), ip' = (p{t,x,v'), p'^ = p{t,x,ii). 

The post-collisional velocities (w,f*) and the precoUisional ones {v',v'^) satisfy 
the conservation of momentum and kinetic energy 

j V + ^ v' + v'^ 

since the collisions are assumed to be elastic. The solutions of (5) can be 
parametrized by a unit vector n G S^^^, so that 

{v' = V — [{v ~ v^).n]n 
v'^, = f* + [{v — v^).n]n. 

The other two collision operators Qi{p) and Q2{p) both take into account the 
collisions between particles from the normal fluid and from the superfluid. The 
operator Qi{p) writes 

Qi{p){v)= I [p{v')p{v:)~p{v){p{v') + ^{v:) + l)]Av'dv: 

JCi 

[fiK) ifiv) + p{v*) + 1) - p{v)p{v^)]dv^dv'^ 

(6) 
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where Ci is the set of all the velocities v' , v'^ such that 

,./l2 , I../ 12 C^) 



+ w = u' + 



and C2 is the set of the velocities , v'^ satisfying 

As to the last collision operator Q2; it writes 

where C3 is the set of all velocities w*, f f ^ such that 

\ -2dte-2Vi + \v,\^ = \v'\^ + K\^. 

The terms —2dt9 — 2Vi which arc added in the microscopic energy conservation 
equations are here to ensure the possibility of having a condensate with velocity 
and energy which are other than zero. Indeed, the usual Boltzmann-Nordheim 
equation was derived in a space homogeneous setting, and up to translations the 
energy and the velocity of the condensate can only be zero in such a framework. 
The potentials Vi and V2 are given by 

f Vi=a{2p„+ps) 

\ V2 = ^{2pn+Ps) 

where a is a positive constant, ps is given by (2) and 

Pn = / (p{t,x,v)dv 



is the mass of the normal fluid. 

Note that very few is known mathematically about the Boltzmann equation 
for the bosons. It has been studied only in the spatially homogeneous case; see 
for example [11], [5], [6]. 

In this work, we will (formally) perform the compressible Euler limit on the 
system (3) together with a semiclassical limit. To do this, we rescale time and 
space by changing {t,x) to (|-, |), where e is a small parameter, which gives the 
scaled system 

dt^ + v.\/^ip - V^yi.V„v? - -0(v3) + —Qi{ip) 
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where h = eh, and we want to determine the asymptotics of this system as 
h ^ and £ ^ 0. 

In Section 2, we detail the main physical properties of the model. Then, in 
Section 3, we obtain a fluid description of the Bose-Einstein condensate and the 
non-condensate phase. We therefore study the fluid system. In Section 4, we get 
some conditions for the system to be hyperbolic, and then we show the existence 
of some local strong solutions. In Section 5, we finally give some details about 
weak solutions and the propagation of shocks. 



2 Properties of the model 

The possibility of obtaining a fluid description of the Bose gas governed by (3) is 
a consequence of some fundamental properties of this system. First, like for the 
classical Boltzmann and Schrodinger equations, we can write the macroscopic 
equations wich govern the mass, momentum and energy for both parts of the 
gas. Then, we characterize the distributions ip which minimize the entropy and 
thus are expected to be attractors of the dynamics in the fast relaxation limit. 
These two informations will allow us to perform the compressible Euler limit in 
Section 3. 



2.1 Conservation laws 

Integrating the kinetic equation in (3) against 1, v, we obtain the equations 
on the mass, momentum and kinetic energy of the normal fluid. From the Gross- 
Pitaevskii equation, we derive the equations for the mass, momentum and energy 
of the condensate. 



2.1.1 Mass 

We plug the expression of (2) in the second equation of the system (3). After 
dividing by the exponential, we can write two equations, one for the real part 
and one for the imaginary part. They read 

dtPs + '^x-iPsUs) = PsQ2{v) (10) 

and 

a*^--^A,V^-^-v^2 (11) 

where 

denotes the bulk velocity of the condensed phase. Equation (10) gives the law 
governing the mass density of the condensate. 

For the kinetic equation, let us first remark that 

V J Qi{^){v)dv = - \ Us \Q2{^)- (12) 
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Then, integrating with respect to the velocity over the whole space M.^ , we 
obtain the equation satisfied by the mass of the normal fluid: 

dtPn + ^x-{PnUn) = -psQ2iv) (13) 

where we defined the bulk velocity of the normal fluid as 

1 



vf{t, X, v)dv. 

Pn Jrw 

We can see that neither p„ nor ps are conserved quantities, unless the thermo- 
dynamic equilibrium is reached, that is, 

Q2iv) = 0. 

However, the total mass 

p = Pn + Ps 

satisfies a conservation law without any assumption of equilibrium: 

dtP + Va:-iPnUn+ PsUs) ^0- (14) 

This means that the normal fluid and the condensate exchange mass until they 
are at thermodynamic equilibrium. 

2.1.2 Momentum 

Integrating the kinetic equation against v, we get 

9t(p„W„) + Vx. I / V^Vipdv] ^-pn'^xVi- PsUsQ2{'f)- (15) 

For the condensed part, we differentiate equation (11) with respect to x, and, 
thanks to the mass equations (10) and (13), we can write 

dt(psUs) + VxXPsUs ® Us) = -psVa;V2 + PsUsQ2{<P) 

Plugging the expression of the potentials (8) , we can see that the total momen- 
tum satisfies 



dt{PnUn+ PsUs) +^ X- \ f V ® Vifidv + PaUg ® Us + a{2pn + Ps)In 
= — Vj;. {^/p^Vly/p^ - Vxy/fh® ^xy/fh) 

where Jjv is the N x N identity matrix. 

The momentum of each fluid is not preserved even at thermodynamic equi- 
librium; however, the total momentum obeys a conservation law. This means 
that, even at the equilibrium, the two parts of the fluid exchange some momen- 
tum. 
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2.1.3 Energy 

We now integrate the kinetic equation against 



For the condensed part, the energy writes /i^jVx^P, and satisfies 



Computing 

and since 

it foUows that 



/i3sR(iV,(A,*).V,*) 



(16) 



+ YV^p,.V:,Q2(<y3) + h'\'^.x^\^Q2{ip) 

On the other hand, the equation for the potential energy writes: 

5t^(2p„ + Ps)^+Va;.(a(2p„ + Ps){2pnUn + PsUs)) 

= a(2p„U„ + PsUs).Vx['2.Pn + Ps) - aps{2pn + Ps)Q2{ip) 
= 2p„U„Vj:Vi + 2psUsV^V2 - PsViQ2{ip) 

thanks to expressions (8) of the potentials. If we add this equation the the 
equation for the energy of the condensed part (16), we note that we obtain 
a conservation equation, plus some terms coming from the interaction with 
the normal fluid. This conservation property is reminiscent from the energy 
conservation of the solutions to the Schrodinger equation. 

Collecting all the pieces together, we get the full energy equation: 



\vYipdv + ps\u,Y + h'\^xVp:\ + 2 + PsY 



v\v\'^ipdv + Ps\Us\'^Us + a{2pn + Ps){2pnUn + PsUs) 



(17) 



h'^T{ip,ps,Us) 
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with 

T{(p,ps,Us) =VJ;.(\/p^(32('^)Vx^^o^) + 5S7x-{'Vx^/pI<E)'^x^/pI : u,^) 

- Vx-{y/fhUsAxy/fh) + '^x-i^x^^'Vx-iy/fhUs)). 

We check that the total energy is conserved, and that all the terms which 
do not have a clear physical interpretation formally vanish in the limit h ^ 0. 

2.2 Thermodynamic equilibrium 

Letting £ ^ in (9) suggests that, in the fast relaxation limit, (p satisfies, for 
almost every v £ M.^ , 

Qiifi) + PsQi{(p) = 0. 
We therefore need to characterize such functions: 

Proposition 1. Let f = ip{v) be a smooth, positive valued function on . Let 
Us G M^, Vi £ M. be some given constants and let 

Pn = [ ip{v)dv, Un = — [ Vip{v)dv. 

Define Q{f) and Qi{Lp) by (4) and (6). Assume that h = 0, so that C2 is the 
set of all the velocities , v'^ such that 

{Us + v'^ = V + 
\us\^ + Vi + \v:\^^\v\^ + \v,\^. 

Then, 

Q (if) + PsQ lif) = 



if and only if ip is a bosonian, i.e. 



1 - M 
where 

M = e 2T . 

for some positive constant T. 
Proof. The identity 

Q{ip) + PsQiiip) = 

implies that 



^(-) = T^, (18) 



{Q{ip) + PsQiiip))\og[j^]dv^O. 
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Using the symmetries of the cohision integrals, we have 



1 + v? 

(1 + ^)(1 + ^,)(1 + ^')(1 + ^'J 



1 + ip' l + ip'^ 1 + <^ 1 + 



and 



(l + ^)(l + V5,)(l + ^'J 



1 + (^t 1 + 1 + 



But for all x, ?/ > we have 

(a;-?/)(logy-loga:) < 

with equality if and only if x = y. Hence both integrands are nonpositive 
functions, and wc claim that: 

Q{ip) + PsQi{ip) = 

if and only if 

Qif) ~ a.e. and Qi{f) = a.e.. 
For the same reason, Q{^p) vanishes if and only if 

l + (p' 1 + l + Lpl + (p^ 

for almost every € K^, n e S^^^, with v',vl satisfying (5), and Qi{(p) 
vanishes only when 

for almost every v € and w*, G C2 with h ~ {ug € M.^ and Vi are fixed). 

We therefore have to solve the functional equation (19). Thanks to the 
classical Boltzmann theory, we know that if for all w, f*, w', satisfying (5) we 
have 

.//* ~ / 

then / is a Maxwellian distribution, namely 

J = — g(a|'Up+fc.«+c) 
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with a, c G R and b E M.^ . Taking / = gives that Q{(p) vanishes if and only 



M 



-+s 



1-M 

The Maxwellian function M can be rewritten 

M = e'- 

for some T > and (5 e M. using the macroscopic velocity, since 

vifdv = p„u„. 

We can compute S by plugging this into Qi{(p) = 0, which leads to 

Ml = MIvU 
for aU V eMJ^ ,v^,v'^ E C2. We thus obtain: 

— = H —=i 



2T 2r 2T 

which gives, thanks to the definition of C2, 

which is the expected result. □ 



Remark 1. M satisfies 

so that 

In particular, 



l + (p 
0<M <1. 

\us - u„p < a{2pn + Ps)- 



Remark 2. In [6], it is shown that the equilibrium solutions of the Boltzmann- 
Nordheim equation, which are the minimizers of the entropy, are the bosonian 
distributions (18) to which a Dirac function must be added if the mass of the 
initial data is greater than some critical value. In fact, for initial data having 
a big mass, the formation of a singularity taking the form of a Dirac function 
means that a condensate has been created in the fluid. Our goal here is to study 
a model (3) where a condensate has already occured, and then the corresponding 
mass has been deleted from the normal component of the fluid. Hence, the initial 
mass of the normal part of the fluid is assumed to be less than the critical value, 
and consequently the equilibrium solution is made only of the regular part. 
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3 Formal derivation of the fluid model 

Our goal in this section is to derive a fluid description of the Bosc gas governed 
by (3). This wiU be achieved at formal level in three steps: we first perform 
the semiclassical limit h 0. It allows to pass from a quantum to a classical 
description of the condensed part. Then, rescaling the time and space as in (9), 
we perform the compressible Euler limit by the moment method. We get a set 
of equations on the masses and momentum of the two fluids, and an equation 
for the total energy. This system is closed. However, it is a system of 2N + 3 
equations and as much unknowns. To simplify the subsequent study, we replace 
it in the third step with a smaller system using the isentropic approximation. 

3.1 Semiclassical limit 

When performed on the usual Schrodinger equation, the semiclassical limit h —^ 
allows to recover the laws of classical mechanics from quantum dynamics [3]. 
Indeed, h can be seen as a measure of how far we are from classical mechanics. 
In our case, the limit will allow us to treat the condensed part as a classical 
fluid. 

We assume that all the quantities converge in suitable functional spaces as 
/i — !■ 0. Moreover, we still write Ps,Us, 9, ip, p„, u„ the limits of these quantities. 
Under these assumptions, the equation for the phase (11) becomes 

Hence, as /i ^ 0, Q2{'f) is changed into 

Q2M = / [^{v'Mv',) - ^{iu){if{v') + + l)]dv,d«'d< 

where the new C3 is the set of all velocities i;*, w', v'^ such that 

Us + = v' + v'* 

\us\^-Vi + \v,\^ = \v'\^ + \v:\^. 

Substituting Q2{^) with its new expression, the mass equations (10), (13) and 

(14) write the same as before. The momentum equation for the normal fluid 

(15) is also unchanged. The equation for the momentum of the condensed part 
becomes: 

dt{psUs) + Vx-{PsUs ® Us) = -ps'^xV2 + PsUsQ2{'P)- 

Thus, the total momentum satisfles the conservation law 

dtiPnUn + PsUs) + V.^,. ( / V® Vifdv + PgUg ® Ug + Q;(2p„ + Ps)In ] = 0. 
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Letting /i ^ in the energy equation (17) leads to: 



2 



+ ^1^1 + Ps\Us\ Us + a{2pn + Ps)(2pnUn + PsUs) j =0. 

Note that in the scmiclassical Umit, the total energy is conserved. 
3.2 Hydrodynamic limit 

We now perform the change of scale {t,x) i-^ {-,-) where e ^ is a small 
parameter. As can be seen on the system (9), this scaling leaves the equations 
associated to the condensed part unchanged, excepted the collision term Q2{^) 
which gets a factor -. The small parameter e is known as the Knudsen number, 
which is proportional to the mean free path, that is, the average length covered 
by the particules of the gas between two collisions. Letting this parameter 
go to means that the collisions occur on a time scale which is very small 
compared with the observation time scale, so that one can consider that the 
local thermodynamic equilibrium is reached almost instantaneously. 

The principle of the moment method for the compressible Euler limit is as 
follows: multiplying the kinetic equation in (9) by e and letting e — ^ suggests 
that Lp converges to the bosonian (18). We thus replace by its limiting value 
in the moment equations, namely, in the equations for the mass, momentum 
and kinetic energy of the normal fluid. We do the same for the equations which 
govern the mass, momentum and energy of the condensate, and we obtain a 
closed set of hydrodynamic equations describing the macroscopic evolution of a 
Bose gas. 

We assume that all the macroscopic quantities converge as e ^ in suitable 
functional spaces. The only technical point is to compute some integrals, in 
particular those associated with the momentum and heat flux. We have the 
identity 

V ® V-^ —dv = PnUn M„ + plN 



defining the pressure 



1 /■ , ,2 /3e 



1 — /3e 2T 



where 



fi = e 2T and < /3 < 1. 



Using obvious symmetry properties, we deduce that 

\v?r^^J^Av = Pn\Un? +Np. 
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In the same way, we have 



M 



1 -M 



dv ^ p„|u„pM„ + (N + 2)pUr. 



At the end, the hydrodynamic limit of the system (3) is given by the following 
set of equations: 



dtPn + ^x-{PnUn) = 
dtPs + Vx-{PsUs) = 

dt{PnU„) + Vx-{PnUn ® +PIn) = -Q!p„V:r(2p„ 
dt{psUs) + Vx-iPsUs ® Us) ^ -|/9sVj;(2p„ + Ps) 
dtaPnWnl^ + ip.>sP + f P + f (2p„ + 



Ps) 



(20) 



■» ' 2 

-f(2p, 



^Ps\Us\'^Us 



iN+2] 



+ Ps){2pnUn + PsUs) 







This is a kind of two-phases Euler system, the second fluid (the superfluid) 
being pressurcless. They do not exchange mass, while they exchange some 
momentum, contrary to what occurs in the model proposed by Landau [10]. 



3.3 The isentropic approximation 

The system (20) is somehow very complex insofar as it is not a system of con- 
servation laws, and defining weak solutions for such a system is complicated. 
Consequently, we approximate it with a simpler system, which can be written 
as a system of conservation laws. To do this, we first claim that the kinetic 
equation in (9) is endowed with a natural entropy S{(p), which is defined by: 



S{^)^—[ ((l + ^)log(l + v^)-(p 

Pn JK" 



— iplogip) dv. 

We will show that in the case when that entropy is approximatively constant, 
we can compute the pressure law, thus reducing the number of unknowns, and 
obtain a closed system of conservation laws. 
The entropy satisfies the following equation: 



dtiPnSm 



w ((1 + (yj) log(l + ip) - ip log (p) dv 



{Q{^) + pMv))\og 
e J^N \l + 'P 



(21) 



dv > 



which is analogous to Boltzmann's H-thcorcm. When e ^ 0, the density func- 
tion (fi goes to its equilibrium value and the entropy (now called S) becomes 



PnS 



M 
1 - M 



logM + log(l - M) 



(22) 
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Then, it comes 



77 ^ TF log "i r? = PnOMn 



^ 1 - Af ° I- M I- M ° 1 - il/ 
and 5 satisfies 

dtiPnS) + Vx-iPnSUn) = 0. 

Moreover, thanks to the mass conservation for the normal fluid (first equation 
in the system (20)), we get 

dtS + Un.VxS = (23) 

which is a transport equation. Note that S = So, where Sq is a constant, is a 
trivial solution of (23). 

Remark 3. In the fast relaxation limit e 0, the entropy dissipation is expected 
to concentrate on shock solutions, as occurs for the classical Boltzmann equation. 
Hence, inequality (21) is expected to become an equality in the limit as long as 
smooth solutions are considered. In other words, (23) should he a combination 
of the other equations of the system. 

Using (22), we can compute the entropy with respect to pn,p,P and T: 
S = iV-i— - log/3 / log(l - (3e- — )dv. 

Pn Pn JRN 

From the expressions of the density and pressure for the normal fluid: 

Pn-T-f^f ^^d., V-^T-'^-'l kP^^d., 
remarking that, thanks to an integration by parts, 



\v\ rrprdf = --T/ log(l-/?e ^ )dv, 

l-/?e-4- P Je" 

we get 



2 J Tpn 

If we write 

p„^T^/^Foil3), p^l-T^/^+^F,i(3), (24) 
we obtain at the end 



2 Nj Fo(/3) 



14 



It comes out that the entropy depends only on the variable (3. If we show that 
/? I— > S{(}) is one-to-one, we will be able to say that if 5' remains constant, so is 
/3; hence, thanks to (24), we will obtain the following pressure law: 



Af + 2 



P = CNPn" ■ (25) 

with 

7V(Fo(/3o))i+2/^ 

for some constant Pq e (0, 1). 

Lemma 3.1. The function f3 G (0, 1) i-^ S{f3) is one-to-one. 
Proof. We compute the derivative of S with respect to /?: 

Hence, 



2 NJ FoiPy (3 



I , 1\ 1 ( f 1^ ,2 e"^ J f f3e' 



e 2 



/3' 



We want to show that S'{(3) is negative. To do this, we expand the integrands 
into power series. We get for the simplest one 

/ ' . dv= V ^ di; = V/3" / e-("+i'^d«, 

that is. 



Wc compute the other integrals using the same tool, and the entropy derivative 
becomes 



■X 



/3 
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We can calculate the products of the summations: 



2 



with 



and 



with 



-E 



(fc + l)^/2(n-/c + l)^/2' 



E („ + l)N/2+l E ^" + l)W/2-l E 
n=0 ^ ' n=0 ^ ' n=0 



s» = E 



A simple calculation shows that 
But since 

it comes 

< 0, 

and thus /? i-^ <S'(/3) is one to one. □ 

From now on, wc will assume that the dimension of the space is iV = 3. 
Moreover, we assume that the gas is translation invariant along the y- and z-axis, 
and that the y- and z-components of the velocities u„ and Us are zero. Under 
these conditions, the system under consideration becomes one-dimensional 

dtPn + dxiPnUn) = 

dtPs + d^iPsUs) = ^25) 

dtiPnUn) + dxipnUn + ^Pn^) = -~apndx{2pn + Ps) 
dtipsUs) + dx{psul) = -jPsdxi^Pn + Ps) 

where c = £3 . 

If {p„, ps,u„, Us) is a regular solution of (26), then ps,Un, Us,p= cpn^) 
will be a regular solution of (20). However, this is not true any more when 
dealing with weak solutions. The reason is that the entropy docs not remain 
constant along shock waves as mentioned in remark 3. Hence, system (26) is 
actually an approximation of system (20), and not just a particular case of it. 
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4 Hyperbolicity and regular solutions 

In view of giving some mathematical justifications to the derivation of the hy- 
drodynamic model (20) and of its isentropic approximation (26), it is necessary 
to have a good understanding of its structure and a good mathematical theory 
of existence and stability. We are therefore interested in this section in finding 
smooth solutions to (26). 

In Section 5, we will deal with weak solutions of this system. As it will be 
more convenient to work with a system of conservation laws, we introduce right 
now the following system of equations: 

dtPn + dxiPnUn) = 

dtPs + dxiPsUs) = 

dtiPnUn + PsUs) + OxiPnU^ + PsUl + Cp^J^ + f (2p„ + Ps^) = 
dti^PnUl + ^Psul + ^Cp^J^ + a{pn + \PsY) 

+dx{^Pnul + \psWl + ^Cp^J^Un + Q;(p„ + \ps){2p„Un + PsUs)) = 

(27) 

Actually, the inversibility of the function 



F 



(Pn, Ps,Un,Us) 



(whose jacobian is jp — pnPs{u 
0+ 



u = 



(Pn\ 
P 

\Us) 



and 



O-^ lU = 



(Pn\ 

P 

\UsJ 



{Pn, PsiPnUn + PsUs, ^PnU^ + ^Ps^s 



3~ I 



a{p„ 
Un)) on the sets 



1 



Psf) 



s.t. PmPs > and w„ > Ug 



s.t. pn,Ps>0 and u„ < 



is enough to say that on these sets, the smooth solutions to (26) are also solutions 
of (27), and conversely. In addition, they are also solutions of 



dtU + Ad^U = 



(28) 



where 



A 



( Un p„ \ 

Us /9f 

cpn ^ + 2a a Un 

\ a % Q uj 



(29) 
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(with c = |c > 0), and any of this problem is (strictly) hyperbolic provided the 
others are. 

4.1 Hyperbolicity 

For a system of the form (28), the right property to study is hyperbolicity, 
that is, the possibility to diagonalize the matrix A given by (29). Indeed, this 
property allows to claim the existence and uniqueness of smooth solutions to 
(28) (see [4]). It will also be useful when dealing with weak solutions. 

Since system (26) is hyperbolic if and only if system (27) or system (28) 
is, it is enough to investigate the eigenvalues of the matrix A given by (29) 
to conclude on the hyperbolicity of any of these systems. The characteristic 
polynomial of A writes 

Even if it is of degree 4 and there are explicit formulae for the roots of polyno- 
mials of degree 4, they are very complicated and we are not able to use them 
here. Our analysis relies only on basic analytical tools to find cases where the 
polynomial has 4 different roots. The main result of this section is the following: 

Theorem 4.1. Systems (26), (27) and (28) are strictly hyperbolic if Pn,Ps > 
and if one of these conditions holds: 

(Un-Usf < Cpl/^ (30) 
a 2/3 

(m„ - Us) < (31) 

cpn + 2ap„ 

Obviously, the condition (32) is the most interesting, because it can be prop- 
agated by the solution, whereas the other two may fail to hold after a small time. 

Proof. We work with U = {pn, Ps^UmUs) and a fixed, with pn,Ps,ce > 0. We 
first notice that 

lim Pa (A) = -foo. 

A — '±00 

On the other hand, if we let 




we can remark that 

,2 



Vi = l,2 P{X*) = -a'p„ps<0. 
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But Us e (AjjAj) satisfies 

Hence, if we assume tliat (30) holds, we fiave 

Pa{us) > 0. 

Since Pa is a continous function of A, we claim that it has four distinct roots 
(Aj;)j-i..4 depending on U and a and satisfying (see Figure 2) 




Figure 2: In the case when Pa{us) > 0, the polynom Pa has four distinct roots. 
In the same way, we let 

{A3 = - cpl!^ + 2apn ^ 
A4 = Un + \l Cprl'^ + 2ap„ 

then 

Pa(\1) = Pa{K) = -a'^PnPs < 
and ^ 

-PA(Wn) = C-Ps/Of/^ " (U„ - Us)^{cpf-/^ + 2apn). 
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Thus, under the condition (31), Pa has four distinct roots (Ai)i=i..4 wich satisfy 



2 / 2/3 

Ai < u„ - y cpn + 2ap„ < A2 < u„ < A3 < u„ + y cpn + 2ap„ < A4. 

We have then proved that under the condition (30) or (31), the systems (26), 
(27) and (28) are strictly hyperbohc. 

Wc now assume that neither (30) nor (31) are verified, that is, 

(w„ - Us)' > cpl/^ (33) 

and 



iun~u,r>^fP^^. (34) 



Thus, Un ^ Us, PA{un) < and Pa{us) < 0. Let us study the sign of Pa( ""^"° )■ 
Pa (^ ""^"' ^ = ^K-"s)*-J + cp2/3 + 2ap„) (^„-u,)2 + c|p,p2/3_ 
Let Q{X) be the polynomial defined by 

Q{X) -^qX'-\ (f + cpy^ + 2«p„) X + c'^PspT- 
Its discriminant is 

^=Y^ (f Ps - cp^/' + 2ap„)' + c|pf/3 > 0, 
and its two roots write 



Hence, 



^- = 3^(i (fp«+CPn/'+2«P„)-V^ 



^ ( " o ' ) > ^ (u„ - Us)' > /?+ or (u„ - usf < (3- 



2 

Let us assume that (32) holds, or equivalently, that 

2apn < 

Then, 

/3+ < ^ (f + 3cp^^ 
But, thanks to (33) and (34), we know that 

(u„ - u,)2 > cp2/3 
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and 

from which we deduce that 



(u„ - Us) > —ps, 



Moreover, we notice that Pj^{us)P^{un) ^ 0. We can conclude that the systems 
(26), (27) and (28) are strictly hyperbolic under the condition (30), or (31), or 
under the only condition (32). In this case, the eigenvalues of A satisfy 

Ai < Un < M < 7: < A3 < lis < A4 



if Un < Us, or the same thing inverting w„ and Ug if Ug < Wn- 



□ 



4.2 Entropy 

We are now interested in finding a convex entropy to system (26), that is, a 
convex function E{U) such that there exists a function G{U) which satisfies 

dtE{U)+d.MU)^0- 

Our research is motivated by the fact that a strictly hyperbolic system having 
a strictly convex entropy admits regular solutions (see [4] for example). 
The entropy wc arc looking for is given by the physical energy: 

EiU) = ^Pnul + ^Psui + Y^Cp?^ + a{p„ + ^Psf ■ 

We can find the equation satisfied by E just from (28); it writes 

dtEiU) + d^GiU) = 

with 



1 



1 



1 



G{U) = 2'^""" + 2*°*"^ ^ 2^PnUn + a{pn + -Ps)(2p„W„ + PsUs) 

Let us now study its convexity. Its Hessian matrix is given by 



cpn + 2a a Un 

a ^ u< 

Un p„ 

Us psj 



If we denote X 



1"^, where is the transposition operator, then 



Us 

X2 

Ps 



Pn 2-3 H Xi 

Pn 



a 

2 " 

CPn ■ 



X2 



2aps 



2a 



aps ~ 2uj 

^ 2aV, 

P- 



■Xl 



aps 



2ui 



Xl- 
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Hence, E will be strictly convex uniformly on every compact subset of O under 
the conditions 

2 Ps 

cpn ^ + 2a - ^ '-^ > 

Pn aps - 2uj 

which are equivalent to 

{2 Q; 
' 2^' . (35) 
Pa{0)>0 

We can now write the following theorem of existence of solutions: 
Theorem 4.2. Let a > 0, and 

A=^^{pn,Ps,Un,Usf, pn,Ps>0, (30) , (3 1) or (32) holds, 

ul<^Ps, Pa(0)>0, Un^U, 

Let Uo E C^{R■,^1) where is a compact subset of A, such that dxUo G 
for some I > ^. Then there exists T^o, such that < T^o < oo, and a unique 
function U G C^{[0,T^) x R;A), which is a classical solution of the problem 
(28) on [0,Too) with the initial condition 

U{Q,x) ^Uq{x), VxeR, 

and of the problem (27) with same initial condition. Moreover, 

dxU{.,t) e c°([o,r^);i/'). 



5 Weak solutions and the Riemann problem 

In Section 4, we have proved the existence of smooth solutions to the problem 
(26), or equivalently to (27) or (28). However, these solutions are not satisfac- 
tory. First, they exist only for small times, while we would like to have global 
solutions. Moreover, physical experiments often show non regular behaviours 
to physical problems. The solutions of hyperbolic systems having a physical 
meaning should therefore be authorized to be discontinous. It is then natural 
to define a weaker notion of solutions. 

In this section, we will recall the usual definition of weak solutions to a 
system of conservation laws, and state an existence and uniqueness result for 
the system (27) . We will then give some details on the qualitative behaviour of 
such a solution by studying the characteristic fields and the shock curves. 



22 



5.1 Existence of weak solutions 

When dealing with weak solutions, systems (26), (27) and (28) are no more 
equivalent. Indeed, the usual definition of weak solutions for an hyperbolic 
system is valid only for systems of conservation laws. We consequently restrict 
our analysis to the system (27), which we rewrite here: 

dtPn + dx{pnUn) = 

dtPs + dxiPsUs) = 

dtiPnUn + PsUs) + dx{pnU^n + Ps^l + CpV^ + f (2p„ + PsY) = 
dti^PnUl + \psUl + \cp^rl^ + Q;(p„ + \ps)'^) 

+dx{^PnUl + \psul + ^CplI^Un + a{pn + ips)(2p„M„ + PsUg)) = 

(36) 

With obvious notations, we can write it 

dtF{U) + d.,H{U) = 0. 
We define weak solutions in the usual way: 

Definition 5.1. Let Uq = (p° , p°, u°, u°) e {L°°{R))^. We say that the bounded 
vector field U = {pn, Ps,Un,Us) is a weak solution of system (36) supplemented 
with the initial condition Uq if for all $ G (C"'^(R-|- x M))** compactly supported 
we have 

f f (F{U{t,x)).dt<i>{t,x) + H{U{t,x)).d^<!>{t,x))dxdt 
i^(f7o(x))$(0,x)da; = 0. 



Moreover, we say that U is a solution of the Riemann problem if the initial 
condition is of the form 



\Uq ifx>{) 



where are two constant states in W^. 

Wc now state the existence and uniqueness result: 

Theorem 5.2. Let U = {pn, Ps^Un^Us) be a constant state such that pmPs > 
and Un "^Us- Let = {p^, p^ju'^jU^) be a bounded function of x having small 
total variation and such that \\U'^ — C/||l=o is sufficiently small; in particular, 
there exists a constant r > such that p^ > r and p^ > r, and such that 
u° — u° > r almost everywhere, or m° — m° > r a.e.. Then the system (36) 
endowed with the intial condition admits a unique weak global entropic (in 
the sense of Liu) solution. 
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Weak solutions of hyperbolic equations are not unique in general. The en- 
tropy criterion of Liu is a way of selecting physical solutions, and weak solu- 
tions are unique in that class. See [1] for details and more general results on 
uniqueness of solutions to hyperbolic systems. General existence theorems for 
hyperbolic systems can be found in [2] or [7]. Let us just mention that this the- 
orem holds true because of the hyperbolicity condition (32) wich is propagated 
by the solutions. 



5.2 Nature of the characteristic fields 

Let A be an eigenvalue of A, that is, a solution of 



(^{X - w„)^ - cp^ - 2apn^ (^{X - Usf - ^Ps^ - Ct^PnPs = 0. (37) 



An associated eigenvector is 

/ 

X = 



{X-Usf-%Ps 



\ 



aps 



\ 



j-{X-uM^^u,f^^y,) 



J 



q(A - Us) 

We can compute the derivatives of A directly from (37) and then we get, letting 



VrrA = 



/ 2^(icpn^^' + 2a)(fPs-(A-^..)2)\ 

jk{cp^J^ + 2apn - (A - u„)2) 
j^{X~u^){^Ps~{X-Usf) 
\^{X - w,,)(cp'/^ + 2ap„ - (A - M„)2)y 



Thus, we get for the associated characteristic field 



WuX.X = 4 



ig.p-/^ + 2«) (fp.-(A-..)^) 

+ ^ps (^cp^ + 2apn - (A - w„)^^ 
- —(A - ii„)^ ((A - iis)2 - ^ps) 

Pn V 2 / 

cpl 2apn - (A - w„) 



It is easy to see that if 



(A - Unf - Cpn - 2apn > 0, 



(38) 
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holds, then the field is genuinely nonlinear. Since 



PA{un - ycp,\ + 2ap,i) = Pa(w„ + y cp,l + 2apn) = -a^PnPs < 0, 

condition (38) always holds for the fields 1 and 4 associated with the extremal 
eigenvalues. Hence these fields are genuinely nonlinear, and they will induce the 
formation of shock or rarefaction waves. Nevertheless, we have not been able to 
establish anything for the fields 2 and 3. 



5.3 Shock curves 



In this part, we are interested in the shock solutions of the Riemann problem 
for system (27). The solutions we are looking for are of the form 



U{t,x) 



U if X < at 
U+ if X > at 



with ~ {p^ , pf , , uf) the left and right constant states, and cr e M the 
shock speed. They must satisfy the Rankine-Hugoniot equations, which write 
for our system 



Pn 
Ps 



PnUn + PsUs 
PnUl + \psUl + ^Cpl!^ + a{pn + \psY 



2' 

Pn Un 
PsUs 

Pnul + Psul + Cpf/^ + f (2p„ 



PnU^n + iPs^l + \cp^J^Un + Q'(p„ + ^Ps){2pnUn + PsUs). 

[X] =x+ -X-. 



2t-s^s ' 2 

Here, the bracket means 
Let 
and 

Ws = Us — a. 
Then, the first two equations of (39) write simply 



(39) 



Wn = u,-,, - a 



[PnWn] = 
[psWs] = 0. 



(40) 



Hence, if we know the left state U = (p„ , p^ , m„ , uf), the shock speed a and 
the right densities p^,p'^, it is easy to compute the right velocities and uf. 
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The third equation in (39) writes 

[pnUi - PnUnCr + Psul - PsU^a + Cp\l^ + ^ (2p„ + Psf'] = 0. 

Using (40) and letting 

M^^ p-w- = pXwX, Ms= p~w- = ptwt, 

and 



it becomes 



_ 1 _ 1 

Pn Ps 

Ml[r^] + M^Ars] + c[pT] + f [(2p„ + Ps?] - 0. (41) 
In the same way, the fourth Rankinc-Hugoniot relation (39) becomes, using (41): 
Ml [tI] + M3 [r2] + 5Af„ [p2/3] + a(2A/„ + M,) [2p„ + p,] = 0. (42) 

Let 



J{Pn,Ps,Un,Us,Cr) 



/ Pn{Un-(r) \ 

Ps{Us - fj) 

Pn{Un - O? + Ps[Us - a? + Cpf/^ + f (2p„ + Ps? 
Pn{Un - O-)^ + Ps{Us - Cr)^ + |p^^^(u„ - Cr) 
V +f (2p„ + P.s) (2p„(Wn - cr) + ^^(it^ - (t))/ 



Equations (40), (41), (42) imply that the Rankine-Hugoniot relations are satis- 
fied if and only if 

J(p+,p+,u+,u+,cr) = J{p-,p-,u-,u-,a). 

Thanks to the implicit functions theorem, for any left state = {p^ , pj jUj) 
and shock velocity ctq such that Z3(p^ J(p,^, pj, u~, uj, do) is invertible 

(take for example p~ = pj = u^^ = I, — ao ~ 0), there exist some right 
states U^{(j) close to for a in a neighbourhood of cto such that 

J{U+{a),a)^J{U-,a), 

that is, system (36) admits some shock solutions, which can be parametrized 
by the velocity cr. Note that the entropy criterion will select just half of each 
shock curve. 

Such shock solutions can be observed with numerical simulations. Examples 
of numerical results using a Lax-Friedrichs scheme (see [9] for example) are 
given by Figure 3 which represents the quantity pn ■ The figure is made of three 
graphics which show the results of the numerical scheme using a spatial grid 
made of 1000, 5000 and 10000 points. The time step is taken so that 

At 
Ax 
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Figure 3: The mass of the normal fluid p. 



As expected, the solutions show 5 different states separated by 4 waves. This 
contrasts with the solutions of the Euler equations which cannot exhibit more 
than 3 waves in general, 2 waves in the isentropic case. The presence of the 
superfluid part induces therefore some changes even in the behaviour of the 
normal part of the fluid compared to a usual (non quantum) fluid. 



Conclusions 

We have derived a two-fluids model for a quantum gas from a system consisting 
of a Boltzmann-like and a Gross-Pitaevskii equations. Our derivation is formal, 
and it seems to be very difficult to obtain a rigourous proof of convergence, at 
least until there is a good existence framework for the quantum Boltzmann-Bose 
equation. However, we have shown that our hydrodynamic model is well-posed, 
which is the first step in a good understanding of the limiting process. 

The limit system (27) that we propose does not seem to exist in the littera- 
ture. It remains to see whether it provides a good description of the evolution 
of a superfluid, before any further mathematical exploration of its properties. 
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